Power series Schrodinger eigenfunctions for a particle on the torus. 
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Abstract 



O 

o 

psj ■ The eigenvalues and a series representation of the eigenfunctions of the Schrodinger 

'. equation for a particle on the surface of a torus are derived. 

<^ \ Introduction. 



The eigenvalues and eigenfunctions of the Schrodinger equation for a particle con- 
strained to a ring are derived in undergraduate quantum mechanics textbooks as examples 
of quantization from periodic boundary conditions [1]. One might expect that finding the 
eigenfunctions \& for the problem of a particle on a toroidal surface with 
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would follow simply from similar considerations, and, perhaps, even foolishly claim this 
to be true while lecturing to one's math methods students. A Google search on "torus 
quantum eigenfunctions" returns many links to esoterica in compactification, chaotic maps 
and plasma dynamics. Here we address a more prosaic textbook type problem: finding 



(V 2 + fc 2 )* = (1) 



the eigenvalues and eigenfunctions (to the best of our knowledge there are no closed form 
solutions) of the Schrodinger equation for a particle on the surface of a torus. 



Solution Method. 



^ | Toroidal coordinates [2] are defined through the relations 

b ; 

c sinha cosd> , n . 

x = (2a) 



cosha — cos/3 

c sinha sin</> 
cosha — cos/3 

c sin/? 



(26) 

(2c) 
cosha — cos/3 

A simple substitution renders Laplace's equation separable in (a,/3,</>), yielding solutions 
in terms of non-integer Legendre functions [3]. This suggests that the eigenvalue problem 
could be solved by examining the limiting case of two toroidal shells or by setting a, the 
coordinate that loosely speaking acts as the "radial" coordinate in the system given by 
equations [2a-2c] to a fixed value. Unfortunately for the former choice, the presence of the 
energy term on the right hand side of Schrodinger's equation precludes the substitution 



that leads to separability. The latter option gives a differential equation that can be solved 
by the methods employed below, but here we choose to use a more intuitive parametric 
representation defined by 

x = (R + a sin.9) cos0 
y = (R + a sin#) sin</> 
z = a cos6>. 



We assume the simplest possible choice for the Hamiltonian, H = 
[4] 
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(3a) 

(36) 
(3c) 

and take for V 2 



V 2 = o" 

A straightforward calculation gives 
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ds 2 = a 2 d9 2 + [R + a sin 9} 2 d<j) 2 , 



so that 



H = - 1 - 
2 
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(4) 



(5) 



(6) 



a 2 d9 2 a[R + asme]d9 [R + asm9] 2 d<p 2 
Making the standard ansatz for the azimuthal eigenfunction x(0) = exp [imcj)], and defining 



a = ^, j3 = 2Ea , gives for equation (6) 
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a cos 9 dip 



m 2 a 2 



d 2 tp 

m 2 ' [l + asin0]d0 [l + asin0] 2 



^ + ^ = 0. 



(7) 



Equation (7) appears somewhat similar to several well known equations from Sturm- 
Liouville theory [5], and as such, suggests a substitution of the form x = sin6*. That this 
sort of substitution is untenable becomes clear when it is noted that 9 runs over [0,27r]; 
x = sm9 can not cover the entire angular range on the torus. Instead, since ip{9) must 
satisfy i/j(9+2iv) = tp(9), we choose to expand the solution in Fourier series with z = exp [i9] 
as 



E 



(8) 



The Hamiltonian given by equation (7) is invariant under 9 — > ir — 9, so the solutions 
of equation (7) can be split into odd and even parity eigenfunctions. For the Fourier 
expansion given by equation (8) this results in 



(-l) n c- n (9) 
r.„ = !-ir +1 c_„. (10) 
solutions; these cases illustrate the general approach 

d dz_ _d_ 
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even parity — > Co = arbitrary, c n 
odd parity — > cq = 0, c n = ( 



We consider first the m 
without the complications that ensue when m j^O. Letting ^ — > §f ^ in equation (7) 
gives after some algebra the three term recursion relation 



2i 

[n(n + 1) - /?]c n+ i + -(/?- n 2 )c n + [/3 - n(n + l)]c n _i = 0. 



(11) 



The n = case of equation (11) gives 



ci = - c (12) 

a 

for the positive parity solutions and 

c-i = ci (13) 

for the negative parity solutions. 

The recursion relation given by equation (11) approaches some fixed constant at large 
n, hence will diverge for arbitrary z. However, unlike that which occurs for standard issue 
two term recursion relations [4], there is no way to force all higher order coefficients to 
zero. To make the series convergent requires more work. 

Consider the positive parity coefficient cn] the numerator of cat will be a polynomial 
of degree N — 1 in with N — 1 roots. The numerator of cat_|_i with be a polynomial of 
degree iV obtained from multiplication of cn by a factor of (f3 — n 2 ) and addition of a lower 
order polynomial in (3. Suppose the lower order cat numerator is factored into its roots as 

cn « 09 - Ajv)(/3 - \%)...(@ - A^" 1 ). (14) 

Then 

c N+1 ~ ((3 - N 2 )(P - X N )(P - \%)...((3 - X%) + F(P) (15) 



(P-x 1 N+1 )(p-x 2 N+1 )...(p-x% +1 ), 



(16) 



so that given the roots of the N th order polynomial they can serve as test roots for the 
(N + 1) . In practice the lower order roots converge very quickly, and, as n increases the 
roots quickly approach the n 2 limit. Once the roots are determined to a given order, they 
can be fed back into the lower order terms and the series safely truncated at order N + 1. 
Co can then be adjusted to whatever boundary value desired (results will be presented in 
the next section). 

For m^O, equation (7) gives a five term recursion relation 
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(l + ^-)(/3 2 -n 2 ) -m 2 o? 



c n + 



^-(n + l) 2 ]+ 2 



C n +1 + 



T [/3 " (n + 2)2] -f 



c n+2 = 0. (16) 



The three lowest coefficients of the positive parity solutions must obey 



Z "1 Z 

[(1 + ^)(3 2 - m 2 a 2 ]c + 2ta((3 - -) Cl + ^- c 2 = 0. 



(17) 



The recursion relation involving the lowest coefficients of the negative parity solutions 
gives = 0, so for both cases there are no restrictions on the two lowest coefficients. The 
full solution for either parity must be some linear combination of two linear independent 
solutions. 

Again consider the positive parity case (the procedure is of course identical for the 
negative parity solutions). To generate two linearly independent solutions, say ipA and 
ijjBi convenient choices are (cq ,c^) = (1,1) and (c B , cf) = (1,-1). Each choice in turn 
will generate series in /?, and the total solution will be some linear combination of the two 
series. To determine eigenvalues and the proper linear combination, insist that at some 
large value of n, say N, 

Ac^ + Bc% = (18a) 
Aci +1 + Bc% +1 = 0. (186) 

The homogenous pair of equations (18a, 18b) have nontrivial solutions only when the de- 
terminant of the c matrix vanishes. The zeros of the determinant are the eigenvalues. The 
ratio of B to A can be determined for each eigenvalue by reinsertion of the eigenvalue into 
equations (18a) and (18b), and A (or B) then fixed by an initial condition. The procedure 
just described will work for m = states also. 

Results. 

In table 1 we show m = eigenvalues E\o for low lying states as determined by the 
power series method in comparison to a Runge-Kutta solution of equation (7) (a = .5 for all 
results). The differential equation is solved by choosing ip(— ^ ) = 1(0) and ip'(—%) = 0(1), 
integrating forward to ^, then integrating from in the lower half plane to — ^p. Either 
ip (for negative parity) or ip' (for positive parity) for both integration paths must vanish at 
Tj-. It is clear that at least for low lying states, convergence to a good number of decimal 
places is obtained. Tables 2 and 3 show some m = 1 and m = 5 eigenvalues. 

Table 4 shows a comparison between power series values and values obtained by solving 
equation (7) numerically. The values indicated in table 4 are the result of summing to n = 6 
which gave at least four digits of accuracy with the differential equation solution, although 
fewer than six is usually sufficient. Finally, table 5 gives a few normalized wave functions. 
i/j22 was listed specifically to show its resemblance to ^io; it is the next to lowest eigenvalue 
of the m = 2 state, but has qualitatively the same character as the m = ground state. 
This is a topic currently under investigation. 

Conclusions. 

We have shown that a Fourier series method can be used to determine eigenvalues 
and eigenfunctions of (at least one version) of the Schrodinger equation for a particle that 
lives on the surface of a torus. 

This work was initially motivated by investigations into two other areas: constrained 
quantum mechanics [6,7,8] and carbon nanotube physics [9,10]. In both areas it proved 
necessary to determine surface wave functions for low-lying states, and to our surprise, 
analytic forms were not to be found. Because of the importance of toroidal geometry in so 
many branches of physics, it is curious that if closed form eigenfunctions exist (in the sense 



that they could be looked up in the standard references), they do not appear in textbooks. 
Should the closed form solutions be known, their publication would certainly provide both 
a useful pedagogical and practical research tool. 
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Table 1 



N=3 
1.134245 
4.242039 



N=5 
1.122415 
4.054351 
9.077862 



N=10 
1.122288 
4.051724 
9.041071 



DE 

1.122286 
4.051722 
9.041070 



E 3 

Low lying m = eigenvalues from roots of cat 

for the differential equation (DE) solution was truncated when continuity of order ~ 10 
in ip(d) or ip'(9) (see text) was obtained. Asterisks indicate a root does not appear at that 
order. 



Fourier coefficient of tp(9). The eigenvalue 
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N 
E 1 

E 3 

Low lying m 



N=3 

.247927 

1.658615 



Table 2 

N=5 
.249375 
1.662639 
4.485872 



N=10 
.249368 
1.663015 
4.476692 



DE 

.249368 

1.663015 

4.476693 



= 1 eigenvalues as roots of the cn coefficient of the Fourier expansion for 



Table 3 

N N=3 N=5 N=10 DE 

E 1 3.732776 3.705405 3.705427 3.705428 

E 2 ******* 8.855027 8.853639 8.853640 

E 3 14.771896 15.289127 15.164616 15.164615 

Low lying m = 5 eigenvalues. The appearance of a larger root before a smaller root is not 
unusual for higher m values. 



Table 4 

Q _7T 7T r\ 7T TT_ 

u 3 4 u 6 2 

ifj FS .908836 .796182 .270886 -.105468 -.479280 

ipDE .908780 .796192 .270874 -.105471 -.479273 

•021 (9) from the Fourier series expansion compared to the differential equation solution. 



Table 5 

^10=. 1853 - .7413 sin# + .0608 cos2# 
^ 22 =.2667 - .7425 sin# + .0052 cos2# 
i/j 31 =.2671 -.5409 sinfl - .4886 cos20 - .2794 sin3fl 

Some representative (un-normalized) wave functions. Only the larger coefficients are 
shown; all others are an order of magnitude less than those given. 



